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MAGNETORESISTIVE RANDOM ACCESS MEMORY SIMULATION 

Joseph J. Nahas 
BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] This invention relates in general to circuit simulations and more particularly to 
simulating magnetoresistive random access memories (MRAMs). 

Description of the Related Art 

[0002] Computer simulation models have been built for MRAM magnetic tunnel 
junctions (MTJ) of MRAM memories. One such example is given in a patent application 
entitled "Method and Apparatus for Simulating a Magnetorsistive Random Access Memory 
(MRAM), having an inventor Joseph J. Nahas, having a filing date of November 22, 2002, 
having a U.S. Patent Application No. of 10/302,203, and having a common assignee, all of 
which is incorporated by reference in its entirety. 

[0003] One type of MRAM MTJ is disclosed in US patent No. 6,545,906, which is 
incorporated by reference in its entirety. The MTJ disclosed in this patent utilizes multiple 
free magnetic layers to achieve a cell that toggles between states when subjected to a 
sequence of magnetic pulses along two directions. 

[0004] What is needed is a simulation model for MRAM cells with multiple free layers. 
BRIEF DESCRIPTION OF THE DRAWINGS 

[0005] The present invention may be better understood, and its numerous objects, 
features, and advantages made apparent to those skilled in the art by referencing the 
accompanying drawings. 

[0006] Figure 1 is a cross-sectional view of one embodiment of an MRAM memory cell 
having multiple free magnetic layers. 

[0007] Figure 2 is a schematic representation of one embodiment of an MRAM memory 
cell. 

[0008] Figure 3 is a cross sectional view of one embodiment of an MTJ. 
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[0009] Figure 4 shows a set of polarization diagrams illustrating one embodiment how the 
polarization of the free layers are affected by magnetic fields applied to a MTJ during a write 
cycle. 

[0010] Figure 5 shows one embodiment of the magnitude of magnetic fields applied to a 
MTJ in the X and Y directions by applying current to write lines during a write cycle. 

[0011] Figure 6 shows plots of one embodiment of measured values of the resistance of 
an MTJ in both the low resistance state and the high resistance state as a function of bias 
voltage across the MTJ. 

[0012] Figure 7 shows one embodiment of a sequence of two magnetic pulses applied to 
an MTJ during a write cycle. 

[0013] Figure 8 shows an example of average conductance values of physical 
measurements from 16 MTJs as a function of bias voltage for both a high resistive state and a 
low resistance state. 

[0014] Figure 9 shows an example of average conductance values of physical 
measurements from 16 MTJs as a function of bias voltage and temperature for the low 
resistance state. 

[0015] Figure 10 shows an example average conductance values of physical 
measurements from 16 MTJs as a function of bias voltage and temperature for the high 
resistance state. 

[0016] Figure 1 1 is an example of a plot of a polynomial coefficient parameters for the 
curves shown in Figure 9 as a function of temperature. 

[0017] Figure 12 sets forth a flow chart of one embodiment of a process for determining 
coefficients of equations that provide for a conductance value of an MTJ as a function of bias 
voltage and operating temperature for both a low resistance state and a high resistance 
according to the present invention. 

[0018] Figures 13-19 set forth one embodiment of a portion of a MTJ computer 
simulation model according to the present invention. 
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[0019] Figure 20 is a block diagram of one embodiment of a computer according to the 
present invention. 

[0020] The use of the same reference symbols in different drawings indicates identical 
items unless otherwise noted. 

DETAILED DESCRIPTION 

[0021] The following sets forth a detailed description of a mode for carrying out the 
invention. The description is intended to be illustrative of the invention and should not be 
taken to be limiting. 

[0022] One embodiment of the present invention relates to a computer simulation model 
for simulating an MRAM memory cell having an MTJ with multiple free magnetic layers 
(free layers). In one embodiment, the simulation model implements a state machine for 
simulation of the switching states of an MTJ. 

[0023] Figure 1 illustrates a cross-sectional view of one embodiment of an MRAM 
memory cell 103 with an MTJ 100 having multiple free magnetic layers (e.g. free layers 305 
and 309 of Figure 3). Memory cell 103 includes a metal line 102 and a metal line 104, 
running parallel to first metal line 102. As illustrated in Figure 1, an insulator layer is present 
between metal line 102 and metal line 104. In one embodiment, metal line 102 and metal line 
104 may be of different materials. In alternate embodiments, metal line 102 and metal line 
104 may be implemented as one line. 

[0024] Memory cell 103 has an MTJ 100 connected to metal line 104 and a metal line 
106, located below MTJ 100. Also below MTJ 100, but electrically isolated from MTJ 100, 
is a metal line 1 14 which runs substantially perpendicular to metal line 102. In the 
embodiment shown, MTJ 100 is located between substantially orthogonal metal lines 102 and 
114. 

[0025] In the embodiment shown, memory cell 103 includes a transistor 130 formed 
within substrate 101. Transistor 130 includes a current electrode 122 and a current electrode 
124, both formed in substrate 101, and a control electrode (gate 123) overlying substrate 101. 
Metal line 106 is coupled to first current electrode 122 of transistor 130 by interconnect 118, 
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interconnect 108, interconnect 1 10, and interconnect 112. Current electrode 124 is coupled to 
an interconnect 115, which, in the current embodiment, is coupled to Vss. In the embodiment 
shown, MTJ 100, metal line 106, interconnects 118, 108, 1 10, 1 12, and 1 15, and metal line 
114 are all formed within an interconnect region 105 located between substrate 101 and metal 
line 104. 

[0026] A current (e.g. Ij , which is shown as going into the page in Figure 1) applied to 
line 1 14 (Axis 1) below MTJ 100 generates a magnetic field (Hi) at MTJ 100 perpendicular 
to the direction of the current according to the "right hand rule". Like wise, a current (I 2 ) 
applied to line 102 above MTJ 100 generates a magnetic field (H 2 ) at MTJ 100 (which at 
MTJ 100 is going into the page in Figure 1). The strength of the magnetic fields at MTJ 100 
are dependent on a number of factors including the spacing between metal lines 102 and 1 14 
and MTJ 100, the width of metal lines 102 and 1 14, the materials in metal lines 102 and 1 14, 
and the structure of MTJ 100. 

[0027] Figure 2 is a schematic representation of memory cell 103. In the embodiment 
shown, metal line 102 is utilized as a write bit line 102, metal line 104 is utilized as a read bit 
line 104, metal line 1 14 is utilized as a write word line 114, and gate 123 of transistor 130 is 
connected to read word line 223. Transistor 130 is coupled between Vss and MTJ 100. MTJ 
100 is represented in Figure 2 as a resistor. 

[0028] In one embodiment, cell 103 is located in an array of MRAM cells, where write bit 
line 102, write word line 1 14, read bit line 104, and read word line 223 extend to other cells 
of the array. 

[0029] Figure 3 is a cross sectional view of one embodiment of MTJ 100. MTJ 100 is 
made up of a number of layers. Top electrode 303 is in electrical contact with metal line 104 
and bottom electrode 321 is in electrical contact with metal line 106. Free magnetic layers 
(free layers) 305 and 309, fixed layer 313, pinned layer 317, and AF pinning layer 319 are 
made of magnetic materials. MTJ 100 includes a tunnel junction dielectric 311 (e.g. made of 
aluminum oxide) between fixed layer 313 and free layer 309. Free layers 305 and 309 and 
coupling layer 307 (e.g. made of ruthium) form a first synthetic antiferromagnet. Fixed layer 
313, coupling layer 315, and pinned layer 317 form a second synthetic antiferromagnet. 
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[0030] As shown by lines 323 and 325, the magnetic field of free layers 305 and 309 can 
be polarized in a number of directions depending upon the current flowing through metal 
lines 102 and 1 14. As shown by arrow 329, the magnetic polarization of pinned layer 317 is 
fixed in one direction. Also shown by arrow 327, the magnetic polarization of fixed layer 313 
is fixed in another direction, opposite the direction of the magnetic polarization of pinned 
layer 317. In other embodiments, MTJs may have other configurations and/or may be made 
of other materials. 

[0031] Referring back to Figure 2, the resistive value of MTJ 100 is indicative of the bit 
state stored in cell 103. To read the bit state stored in cell 103, a voltage is applied to read bit 
line 104 and read word line 2223. The resistive value of MTJ is then measured by a sense 
amplifier (not shown). The resistive value of MTJ 100 is dependent upon the direction of the 
magnetic polarization of free layers 305 and 309 of MTJ 100. 

[0032] Figure 4 shows a set of polarization diagrams illustrating how the polarization of 
the free layers 305 and 309 are affected by magnetic fields applied to MTJ 100 during a write 
cycle of MTJ 100. Magnetic fields are applied to MTJ 100 by applying voltages to metal 
lines 102 and 1 14 which generate currents in those lines. The magnetic fields rotate the 
polarity of layers 305 and 309 to write a bit value in MTJ 100. Metal lines 1 14 and 102 are in 
one embodiment, simple resistive metal lines with simple linear current-voltage 
characteristics as shown below: 



Ii=G,Vi (1) 
and 

h = G 2 V 2 (2) 

where Ii is the current through line 1 14, Gi is the conductance of line 1 14 in cell 103, and Vi 
is the voltage across line 1 14 in cell 103. I 2 is the current through line 102, G 2 is the 
conductance of line 102 in cell 103, and V 2 is the voltage across line 102 in cell 103. 

[0033] In some embodiments, there is a linear relationship between the current in a write 
line (e.g. 102 and 1 14) and the magnetic field generated at the MTJ by that current. The 
relationships can be written as: 
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H, = K,I, 



(3) 



and 



H2 — K2I2 



(4) 



where Hi is the magnetic field at MTJ 100 generated by a current flowing through line 114, 
Ki is the relationship between the current in line 114 and the magnetic field the current 
generates, H 2 is the magnetic field at MTJ 100 generated by a current flowing through line 
102, and K 2 is the relationship between the current in line 102 and the magnetic field that 
current generates. In one embodiments, the write currents through lines 102 and 1 14 are on 
the order of 10 to 20 mA and the K's are on the order of 5 to 10 Oe/mA. 

[0034] As shown in Figure 4, polarization diagrams 401 405, 407, 409, and 411 each 
include 3 arrows to show the magnetic polarization of fixed layer 313, the magnetic 
polarization of free layer 305, and the magnetic polarization of free layer 309 during a write 
cycle. 

[0035] Figure 5 shows the magnitude of magnetic fields applied to MTJ 100 in the X and 
& Y directions by applying current to lines 102 and 1 14 during a write cycle. Applying 
current through one of lines 102 and 114 generates a magnetic field (H y ) in the Y direction 
and applying current through the other of lines 102 and 114 generations a magnetic field (H x ) 
in the X direction. The X and Y directions are in the same plane as Hi and H 2 in Figure 2. 
The timing of H y and H x during one example of a write cycle as illustrated in the polarization 
diagrams of Figure 4 is shown in Figure 5. 

[0036] Referring back to Figure 4, polarization diagram 401 shows the polarizations of 
fixed layer 313, free layer 305, and free layer 309 when MTJ 100 is in a low resistance state. 
At the start of a write cycle of cell 103, the polarization of free layer 309 is parallel to the 
polarization of fixed layer 313 while the polarization of free layer 305 is in the opposite 
direction. When the Y direction magnetic field (H y ) is applied, the polarization of free layers 
305 and 309 rotate so that their net magnetic field (as designated by arrow 413) in 
polarization diagram 405 aligns with the polarization of the external field H y (as designated 
by arrow 415). When the X direction magnetic field ( H x ) is applied at the same time as H y 
(as shown in polarization diagram 407), the polarization of free layers 305 and 309 again 
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rotate so that their net field aligns with the net external field (see arrow 417). When magnetic 
field H y terminates, the polarization of free layers 305 and 309 again rotate so their net field 
aligns with the external field (as designated by arrow 419 in polarization diagram 409). 
Finally when field H x terminates, the polarization of free layers 305 and 309 again rotate to 
the nearest rest state, which is opposite of their initial rest state, i.e. free layer 309 is polarized 
opposite to that of the fixed layer 305 and free layer 305 is polarized in the same direction as 
the polarization of fixed layer 3 13, as shown in polarization diagram 411. In polarization 
diagram 411, MTJ 100 is in a high resistance state. 

[0037] While a particular write cycle pulse sequence with overlapping magnetic field 
pulses with the Y pulse first, is shown in Figures 4 and 5, the MTJ 100 will respond to other 
pulse sequences. For example, sequential overlapping magnetic field pulses with the X pulse 
first will also toggle the bit by rotating the free layers in the opposite direction to change the 
bit state of the memory cell. However, isolated X or Y magnetic pulses will not change the 
bit state of the memory cell. On the other hand, some pulse sequences such as simultaneous 
or almost simultaneous X and Y magnetic pulses can leave the bit in an unknown state. 
Embodiments of a simulation model are designed to take these various possible state 
transitions into account during a simulation of an MRAM cell. 

[0038] In some embodiments, when the two magnetic layers (e.g. fixed layer 313 and free 
layer 309) on either side of the tunnel junction dielectric (e.g. 31 1) are polarized in the same 
direction, the electron spins states in the materials of both those magnetic layers match, and 
the MTJ is in a low resistance state. However, when those two magnetic layers are polarized 
in opposite directions, the electron spin states do not match and the MTJ is in a high 
resistance state. 

[0039] Figure 6 shows plots of measured values of the resistance of an MTJ of one 
embodiment in both the low resistance state (Ri 0 ), and the high resistance state (R h i) as a 
function of bias voltage across the MTJ. Usually, the resistance characteristics of an MTJ are 
expressed in terms of its resistance in the low resistance state at zero bias voltage (Ri 0 o), its 
resistance in the high resistance state at zero bias voltage (RhioX ar *d the Magnetoresistive 
Ratio (MR). The MR is the ratio of the resistance in the high resistance at zero bias voltage 
to the resistance at the low resistance state at zero bias voltage expressed as a percentage 
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above unity. For the embodiment of an MTJ whose measured results are shown in Figure 6, 
Rioo is approximately 17 KOhms, R h io is approximately 22 KOhms, and the MR is 
approximately 29% as shown in Equation 5 below: 



R io0 nm 



[0040] Since an MTJ includes two metal layers separated by a dielectric layer (e.g. tunnel 
junction dielectric 311), it also has a capacitive component. The capacitance of the MTJ can 
be approximated by: 

A 

C MTJ - E r E o~ ( 6 ) 

where £r is the relative dielectric constant of the dielectric, (AI2O3 in one embodiment), 80 is 
the permittivity of free space, A M tj is the area of the MTJ, and t ox is the thickness of the 
tunnel junction dielectric in the MTJ. 

[0041] Figure 7 shows a sequence of two magnetic pulses applied to an MTJ during a 
write cycle. These pulses are generated by applying current pulses through the write lines 
(e.g. lines 1 14, and 102). For one embodiment of a valid write cycle, the magnetic field of 
the first pulse 703 must reach a higher threshold (H S wi) than a threshold (H S w2) required from 
the magnetic field of the second pulse 705. The higher threshold is indicative of the higher 
rotational force required to rotate the polarization of the free layers from their rest position 
(e.g. as shown in polarization diagram 401) than when they are partially rotated (e.g. as 
shown in polarization diagrams 405, 407, and 409) . In some embodiments, the higher 
magnetic field threshold is required for the first pulse, regardless of whether the first pulse is 
generated by current through line 102 or line 1 14. 

[0042] In some embodiments, the polarizations of the free layers take time to stabilize 
while being rotated. Accordingly, minimum time requirements taul, tau2, and tau3 are 
required in the model between pulse edges as shown in Figure 7. Thus, in some 
embodiments, free layer polarization rotations with shorter times are not predictable and 
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consequently cause errors in the model. Also in some embodiments, an additional minimum 
time (tauO) between write cycle pulse sequences is required to allow the polarizations of the 
free layers to stabilize in their rest state. In one embodiment, the minimum times are on the 
order of 1 to 2 ns. 

[0043] In one example, first pulse 703 is generated by line 1 14 and second pulse 705 is 
generated by line 102 causing the polarizations of the free layers to be rotated one way. In 
another example, first pulse 703 is generated by line 102 and second pulse 705 is generated 
by line 114 causing the polarizations of the free layers to be rotated in an opposite way. 

[0044] In some embodiments of a computer simulation model of an MRAM cell, a state 
machine can be implemented to model a write cycle of the MRAM cell. In some 
embodiments, the states of the state machine not only account for magnetic pulses during a 
write cycle of the MTJ, but also pulses acting on the MTJ generated during a write cycle of 
another MTJ located in the same row or column as the MTJ. In those instances, the MTJ 
being modeled may only be acted upon by a single magnetic pulse. 

[0045] In one specific embodiment of a computer simulation model of an MRAM cell, a 
state machine includes four state variables. The first state variable is designated as "pulse 1." 
In the embodiment of Figure 1, pulse 1 is set (to 1) when the magnetic field generated by the 
lower write line (e.g. line 1 14) passes the higher threshold (e.g. H S wi) if it is the only 
magnetic pulse or the lower threshold (e.g. H S w2) if the other pulse (the pulse generated by 
the upper write line, e.g. line 102) has already started. Pulse 1 is reset (to 0) when the pulse 
goes below the lower threshold (H S w2). 

[0046] The second state variable is designated as "pulse 2." Pulse 2 is set (to 1) when the 
magnetic field generated by the upper write line (e.g. line 102) passes the higher threshold 
(e.g. Hswi) if it is the only magnetic pulse or the lower threshold (e.g. Hsv/2) if the other pulse 
(the pulse generated by the lower write line, e.g. line 114) has already started. Pulse 2 is reset 
(to 0) when the pulse goes below the lower threshold (H S w2). 

[0047] The third state variable is designated as "1-2 pulse train." The 1-2 pulse train state 
variable is used to indicate when a pulse generated by the lower write line (e.g. line 1 14) is • 
the first pulse (e.g. pulse 703 in Figure 7) and a pulse generated by the upper write line (e.g. 
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line 102) is the second pulse (e.g. pulse 705 in Figure 7) in a write cycle pulse sequence. The 

1- 2 pulse train is set (to 1) when a write cycle pulse sequence starts with a pulse generated by 
the lower write line (e.g. line 1 14) (when its the first pulse in a write cycle pulse sequence) 
and a subsequent pulse generated by the upper write line (e.g. line 102) (the second pulse in 
the write cycle pulse sequence) passes the lower threshold (e.g. H S w2). The 1-2 pulse train is 
reset (to zero) when the pulse generated by the upper write line (e.g. line 102) (the second 
pulse in the write cycle pulse sequence) goes below the lower threshold (H S w2), indicating the 
end of a write cycle pulse sequence. 

[0048] The fourth state variable is designated as "2-1 pulse train." The 2-1 pulse train 
state variable is used to indicate when a pulse generated by the upper write line (e.g. line 102) 
is the first pulse (e.g. pulse 703 in Figure 7) and a pulse generated by the lower write line (e.g. 
line 1 14) is the second pulse (e.g. pulse 705 in Figure 7) in a write cycle pulse sequence. The 

2- 1 pulse train is set (to 1) when a write cycle pulse sequence starts with a pulse generated by 
the upper write line (e.g. line 102) (when its the first pulse in a write cycle pulse sequence) 
and a subsequent pulse generated by the lower write line (e.g. line 1 14) (the second pulse in 
the write cycle pulse sequence) passes the lower threshold (e.g. H S w2). The 2-1 pulse train is 
reset (to zero) when the pulse generated by the lower write line (e.g. line 1 14) (the second 
pulse in the write cycle pulse sequence) goes below the lower threshold (H S w2), indicating the 
end of a write cycle pulse sequence. 

[0049] For the specific embodiment described herein, the states of the state variables are 
changed in response to the magnetic pulses generated by the lower write line (e.g. line 1 14), 
and the upper write line (e.g. line 102). 

[0050] When a magnetic pulse generated by the lower write line (e.g. line 1 14) is detected 
as passing the lower threshold ( H S w2) while increasing in magnitude, the results of the action 
are dependent on the states of the pulse 1, pulse 2, 1-2 pulse train, and 2-1 pulse train state 
variables. The Karnaugh Map in Table 1 indicates the results of the detection of this action. 
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Table 1 



Pulse generated By the Lower Write Line (e.g. Line 114) 

Passing Lower Threshold (H S w2> In Increasing Direction 
Where I indicates no action, E indicates an error, and S indicates a state change. 



pulse 1 /pulse2 


1-2 pulse train / 2-1 pulse train 


00 


0 1 


1 1 


1 0 


00 


I 


E 


E 


E 


0 1 


S 

check taul timer 

pulse 1 set 
2-1 pulse train set 
tau2 timer reset 


E 


E 


S 

pulse 1 set 


1 1 


E 


E 


E 


E 


1 0 


E 


E 


E 


E 



[0051] Most of the cells in the Karnaugh map of Table 1 contain "E" (for error) for states 
that should not be logically reachable. For example, any state where 1-2 pulse train and 2-1 
pulse train are both logical Is is not reachable. The "I" (for ignore) indicates that although 
the transition is valid, there is no state change. The "S" indicates a state transition with the 
changes in the state variables provided in Table 1 . 



[0052] When a magnetic pulse generated by the lower write line (e.g. line 1 14) passes the 
lower threshold ( H S w2) with none of the state variables set (state 0,0,0,0), the pulse transition 
is ignored since this might be a "first pulse" of a write cycle pulse sequence. In such case, the 
first pulse must therefore pass the higher threshold (H S wi) to be a valid first pulse of a write 
cycle pulse sequence. 

[0053] When pulse 1 is not set, pulse 2 is set, and neither pulse train is active (state 
0,1,0,0), the transition of a pulse generated by the lower write line (e.g. line 1 14) through the 
lower threshold ( H S w2) indicates the start of a 2-1 pulse train so both pulse 1 and 2-1 pulse 
train variables are set after checking that the taul time requirement in the first phase of the 
pulse train is met. Also, the tau2 timer is started. In one embodiment, the taul time 
requirement is determined to be met if the occurrence of a pulse generated by the lower write 
line (e.g. line 1 14) passing through the lower threshold ( H S w2) occurs at least taul time after 
the pulse generated by the upper write line (e.g. line 102) passes through the higher threshold 
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(H S wi). See Figure 7 for an illustration of the taul, tau2, tau3, and tauO times during a write 
cycle pulse sequence. 

[0054] When pulse 1 is not set, pulse 2 is set, and 1-2 pulse train is set (state 0,1, 1,0), the 
transition of a pulse generated by the lower write line (e.g. line 1 14) through the lower 
threshold ( Hsw2) indicates that there was "ringing" in the lower write line (e.g. line 114) 
current near the lower threshold ( Hsw2) as the pulse generated by the lower write line (e.g. 
line 114) is turning off. Thus, the pulse 1 variable is set. 

[0055] All other states of Table 1 are not logically reachable and therefore are errors. 

[0056] Table 2 is the Karnaugh Map for a pulse generated by the lower write line (e.g. 
line 1 14) passing through the higher threshold ( H S wi) while increasing in magnitude. When 
all the state variable are not set (state 0,0,0,0) a pulse generated by the lower write line (e.g. 
line 1 14) passing through the higher threshold (H S wi) indicates the start of a first pulse of a 
write cycle pulse sequence or single pulse, so pulse 1 is set after the tauO timer is checked. 
The taul timer is started. 

[0057] In Table 2, all other states where pulse 1 is not set are erroneous (indicated by E). 
If in Table 2, pulse 1 is set, then the states are either erroneous or can be ignored since they 
are due to a pulse "ringing" near the upper threshold and do not change the state. Table 2 is 
set forth below. 

Table 2 



Pulse Generated by the Lower Write Line (e.g. Line 114) 
Passing Upper Threshold (H S wi) in Increasing Direction. 
Where I indicates no action, E indicates an error, and S indicates a state change. 



pulse 1 / pulse2 


1-2 pulse train / 2-1 pulse train 


00 


0 1 


1 1 


1 0 


00 


S 

check tauO timer 

pulse 1 set 
taul timer reset 


E 


E 


E 


i 0 1 


E 


E 


E 


E 


1 1 


E 


I 


E 


I 


1 0 


I 


I 


E 


E 
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[0058] Table 3 is similar to Table 1 and Table 4 is similar to Table 3 except that Tables 3 
and 4 apply to a pulse generated by the upper write line (e.g. line 102), the pulse 1 and pulse 2 
variables are interchanged, and the 1-2 pulse train and 2-1 pulse train variable are 
interchanged. 

Table 3 

Pulse Generated By The Upper Write Line (e.g. Line 102) Passing Lower Threshold 

(Hsw2) In Increasing Direction. 
Where I indicates no action, E indicates an error, and S indicates a state change. 



pulse 1 /pulse2 


1-2 pulse train / 2-1 pulse train 


00 


0 1 


1 1 


1 0 


00 


I 


E 


E 


E 


0 1 


E 


E 


E 


E 


1 1 


E 


E 


E 


E 


1 0 


S 

check taul timer 

pulse 2 set 
1-2 pulse train set 
tau2 timer reset 


S 

pulse 2 set 


E 


E 



Table 4 

Pulse Generated By the Upper Write Line (e.g. Line 102) 
Passing Upper Threshold (H S wi) While Increasing. 
Where I indicates no action, E indicates an error, and S indicates a state change. 



pulse 1 /pulse2 


1-2 pulse train / 2-1 pulse train 


00 


0 1 


1 1 


1 0 


00 


S 

check tauO timer 

pulse 2 set 
taul timer reset 


E 


E 


E 


0 1 


I 


E 


E 


I 


1 1 


E 


I 


E 


I 


1 0 


E 


E 


E 


E 



[0059] Table 5 is the Karnaugh Map for a magnetic pulse generated by the lower write 
line (e.g. line 1 14) passing the lower threshold ( Hsw2) while decreasing in magnitude. If all 
the state variable are not set (state 0,0,0,0), the pulse transition is ignored. This would 
normally be indicative of a stand alone, lower magnitude pulse that did not have a high 
enough magnitude to set the pulse 1 variable. 
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[0060] If pulse 1 is set and the other state variable are not set (state 1,0,0,0), a magnetic 
pulse generated by the lower write line (e.g. line 114) passing the lower threshold ( H S w2) 
while decreasing in magnitude indicates the end of a stand alone pulse. In this state, pulse 1 
is reset and the tauO timer is reset. 

[0061] If both pulse 1 and pulse 2 are set and 2-1 pulse train is set (state 1,1,0,1), a 
magnetic pulse generated by the lower write line (e.g. line 1 14) passing the lower threshold 
(H S w2) while decreasing in magnitude is due to "ringing" at the start of the second pulse in a 
write cycle pulse sequence. To maintain the state variables, both pulse 1 and 2-1 pulse train 
are reset. 

[0062] If both pulse 1 and pulse 2 are set and 1-2 pulse train is set (state 1,1,1,0), a 
magnetic pulse generated by the lower write line (e.g. line 114) passing the lower threshold 
(Hsw2) while decreasing in magnitude is due to the end of the first pulse in the write cycle 
pulse sequence. After the tau2 timer is checked, pulse 1 is reset and the tau3 time is reset. 

[0063] If pulse 1 is set, pulse 2 is not set, and 2-1 pulse train is set (state 1,0,0,1), then a 
magnetic pulse generated by the lower write line (e.g. line 1 14) passing the lower threshold ( 
Hsw2) while decreasing in magnitude is due the end of a 2-1 write cycle pulse sequence. 
After a check of the tau3 timer, pulse 1 and 2-1 pulse train are reset, the tauO timer is reset, 
and the bit state is toggled. All other states of Table 5 are not logically reachable and 
therefore are errors. 
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Table 5 

Pulse generated by the Lower Write Line (e.g. line 114) 
Passing Lower Threshold (H S w2) in Decreasing Direction. 
Where I indicates no action, E indicates an error, and S indicates a state change. 



pulse 1 /pulse2 


1-2 pulse train / 2-1 pulse train 


00 


0 1 


1 1 


1 0 


00 


I 


E 


E 


E 


0 1 


E 


E 


E 


E 


1 1 

1 1 ! 


E 


S 

pulse 1 reset 
2-1 pulse train reset 


E 


S 

check tau2 timer 

pulse 1 reset 
tau3 timer reset 


1 0 


S 

pulse 1 reset 
tauO timer reset 


S 

check tau3 timer 
pulse 1 reset 
2- 1 pulse train reset 
tauO timer reset 
toggle bit state 


E 


E 



[0064] Table 6 is similar to Table 5 except that Table 6 applies to a pulse generated by the 
upper write line (e.g. line 102) and the pulse 1 and pulse 2 variables are interchanged and the 
1-2 pulse train and 2-1 pulse train variables are interchanged. 



Table 6 

Pulse Generated by the Upper Write Line (e.g. Line 102) 
Passing Lower Threshold (HSW2) in Decreasing Direction. 
Where I indicates no action, E indicates an error, and S indicates a state change. 



pulse 1 /pulse2 


1-2 pulse train / 2-1 pulse train 


00 


0 1 


1 1 


1 0 


00 


I 


E 


E 


E 


0 1 


S 

pulse 2 reset 
tauO timer reset 


E 


E 


S 

check tau3 timer 
pulse 2 reset 
1-2 pulse train 

reset 
tauO timer reset 
toggle bit state 


1 1 


E 


S 

check tau2 timer 

pulse 2 reset 
tau3 timer reset 


E 


S 

pulse 2 reset 
1-2 pulse train 
reset 


1 0 


E 


E 


E i 


E 
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[0065] In one embodiment, a model for simulating an MRAM cell includes a model for 
MTJ (e.g. 100) conductance characteristics that is dependent upon a voltage across the MTJ 
(bias voltage) and MTJ operating temperature. In one embodiment, parameter values for the 
model are calculated using an empirical model from measured data. 

[0066] Figure 12 sets forth a flow chart for a process for determining coefficients of 
equations that provide for a conductance value of an MTJ as a function of bias voltage and 
operating temperature for both a low resistance state and a high resistance. These equations 
and coefficients may be used in a computer simulation model of an MRAM cell including e.g. 
a model implementing a state machine (e.g. the state machine described above with respect to 
Tables 1-6). 

[0067] In 1201, physical measurements of conductance values are obtained from a 
number MTJ (e.g. 16) over a range of bias voltage increments and operating temperature 
increments for both the high resistive state and the low resistive state. In one example, 
measurements are made for both positive bias voltages and negative bias voltages. 

[0068] Figure 8 shows an example of average conductance values (in mmhos) of physical 
measurements from 16 MTJs as a function of bias voltage for both a high resistive state 
(G R hiu P and G Rh idn) and a low resistance state (Gri g ) at a temperature of 30 C. Regarding the 
high resistance state conductance, the subscript "up" is used for a negative bias voltage 
indicating that the current in the MTJ is "up" through the device and the subscript "dn" is 
used for the positive bias voltage indicating that the current in the MTJ is "down" through the 
device. Also shown in Figure 8 are curve fits (803, 805, and 807) for the averaged measured 
data. For the embodiment, shown, the conductance measurements were made at bias voltage 
increments of approximately 0.10 Volts. 

[0069] As shown in Figure 8,* a single second order polynomial (as shown by line 807) 
provides a very good fit to the data (Gri 0 ) of the low resistance state while two second order 
polynomials, one for positive bias (line 805) and one for negative bias (line 803), are utilized 
for the high resistance state to have a similar level of fit. In one embodiment, the curve for 
the low resistance state is given below in equation 7: 

G = 17.99V 2 + 3.28V + 59.25 (7) 
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where V is the bias voltage in volts and G is the conductance of the MTJ in mmhos. The 
curve for the high resistance state positive bias voltage and the curve for the high resistance 
state negative bias voltage are given below in equations 8 and 9, respectively: 

G = 14.68V 2 + 13.84V + 45.83 (8) 

G = 14.69V 2 - 8.29V + 45.66 (9) 

[0070] In 1203, the conductance measurements for the low resistive state, the high 
resistive state positive bias voltages, and the high resistive state negative bias voltages for 
various operating temperatures are fit into second order polynomial curves. In some 
embodiments, the fittings are performed using the least squares error measure or weighted 
least squares error measure. 

[0071] Figure 9 shows the average conductance values (in millimhos) of physical 
measurements from the same 16 MTJs as a function of bias voltage and temperature for the 
low resistance state (high conductance state). Each curve (913,911, 909, 907, 905, and 903) 
represents a fit of a second order polynomials for each operating temperature at which 
measurements were made. In the embodiment shown, the measurements were made for 
operating temperatures of 30 degrees, 40 degrees, 60 degrees, 80 degrees, 100 degrees, and 
120 degrees Celsius. The curve fits for the measurements at each temperature are given 
below: 



30°C: G = 18.0V 2 + 3.3V + 59.3 


(10) 


40°C: G = 18.3V 2 + 3.4V + 59.3 


(11) 


60°C: G = 18.8V 2 + 3.4V + 60.2 


(12) 


80°C: G = 19.4V 2 + 3.6V + 60.9 


(13) 


100°C: G= 19.6V 2 + 3.6V + 61.2 


(14) 


120°C: G = 21.2V 2 + 3.7V + 62.4 


(15) 



[0072] Figure 10 shows the average conductance values (in millimhos) of physical 
measurements from the same 16 MTJs as a function of bias voltage and temperature for the 
high resistance state (low conductance state). Second order polynomials curves are fitted to 
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the average measurements at each operating temperature for a positive bias voltage (curves 
1029, 1027, 1025, 1023, 1021, and 1019) and for a negative bias voltage (curves 1013, 1011, 
1009, 1007, 1005, and 1003). The curve fits for the measurements at each temperature for the 
- positive bias voltage are given by the equations below: 



30°C: G = 12.6V 2 + 15.9V + 45.5 


(16) 


40°C: G = 13.0V 2 + 15.6V + 46.0 


(17) 


60°C: G= 13.7V 2 + 15.3V + 46.9 


(18) 


80°C: G = 14.7V 2 + 14.8V + 47.9 


(19) 


100°C: G = 15.0V 2 + 14.8V + 48.4 


(20) 


120°C: G = 16.9V 2 + 13.8V + 50.5 


(21) 


[0073] The curve fits for the measurements at each temperature for the positive bias 
voltage are given by the equations below: 


30°C: G = 14.7V 2 - 8.3V + 45.7 


(22) 


40°C: G = 14.9V 2 - 8.3V + 46.0 


(23) 


60°C: G = 15.3V 2 - 8.2V + 46.8 


(24) 


80°C: G = 16.2V 2 - 7.6V + 47.8 


(25) 


100°C: G= 16.4V 2 -7.7V + 48.2 


(26) 


120°C: G= 18.1V 2 -6.7V + 50.3 


(27) 



[0074] Referring back to Figure 12, in 1205, the individual polynomial coefficients for 
the low resistance state conductance curves, the individual polynomial coefficients from the 
high resistance state positive bias voltage conductance curves, and the individual polynomial 
coefficients from the high resistance state negative bias voltage conductance curves are each 
fit to a first order temperature polynomial curve. The conductance curves for each 
temperature are of the form: 

G = A + BV + CV 2 (28) 
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[0075] In 1205, the polynomial coefficient parameters (A, B, and C) for the low 
resistance state, the high resistance state positive bias voltage, and the high resistance state 
negative bias voltage are each made into a function of temperature to reduce the number of 
parameters of a model. 

[0076] Figure 1 1 is a plot of a polynomial coefficient parameters (A,B, C) for the curves 
shown in Figure 9 (equations 16-21) as a function of temperature. As shown in Figure 1 1 for 
the case of the low resistance state, a linear temperature function provides a reasonable fit for 
the polynomial coefficient parameters. The curves for A, B, and C for the data plotted in 
Figure 9 are given below: 



A = .032T + 58.3 (29) 

B = 0.004T + 3.2 (30) 

C = 0.032T+ 16.9 (31) 
where T is temperature. 

[0077] Consequently, substituting equations 29, 30, and 31 into equation 28, the 
conductance curves shown in Figure 9 an be modeled as: 

G = (0.32T + 58.3) + (0.004 T + 3.2)V + (0.032T + 16.9) V 2 (32) 

[0078] A linear fit also applies to polynomial coefficient parameters of the curves of the 
conductance of the MTJ in the high resistance state for both for positive bias voltage (curves 
1019, 1021, 1023, 1025, 1027, and 1029) and for the negative bias voltage (curves 1003, 
1005, 1007, 1009, 1011, and 1013). 

[0079] Thus, the three temperature dependent second order polynomials in voltage can 
each have the form of: 

G = G 0 [(l + C T AT) + (Cvi + CjviAT)V + (C v2 + Ctv 2 AT)V 2 ] (33) 

where Go is the zero bias conductance at the reference temperature, Ct is a temperature 
coefficient parameter of the 0 th order voltage coefficient (A), C v i is the 0 th order temperature 
coefficient parameter for the 1 st order voltage coefficient (B), AT is the change in temperature 
from a reference temperature, Ctvi is a 1 st order temperature coefficient parameter of the 1 st 
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order voltage coefficient (B), Cv2 is the 0 order temperature coefficient parameter for the 2 
order voltage coefficient (C), Cjvi is a 1 st order temperature coefficient of the 2 nd order 
voltage coefficient (C). Using the format in Equation 33 , the conductance function for the 
low resistance state can then be written as: 

G Rlo = GwooKl + C T GrloAT) + (CyiGrlo + GTVlGrioAT)V + (C V2 Grlo + Gxv2GrloAT)V 2 ] (34) 

where G Rlo0 =-j- 

K lo 

where Ri 0 o is the resistance for the zero bias voltage, reference (e.g. room) temperature, low 
resistance state resistance. For equation 34, the subscript "GRlo" is added to the coefficient 
named above with respect to equation 33 to indicate that the coefficients are for the low 
resistive state. 

[0080] The two parts (the negative bias voltage part and the positive bias voltage part) of 
the conductance in the high resistance state must be continuous at zero bias voltage, 
restricting the two polynomials to the same value of G 0 . The two high resistance state 
conductance functions can be written as follows: 



GRhiup = GRhio[(l + CxGrhiAT) + (CviGrhiup + CxviGrhiupAT)V (35) 

+ (Cv2Grhiup + CrV2GrhiupAT)V ] 

G Rh idn = G Rh io[(l + GrGrhiAT) + (CviGrhidn + C T VlGrhidnAT)V (36) 

+ C V2 Grhidn + CxV2GrhidnAT)V ] 

where G Rhi0 = 1 1 



[0081] For equation 35, the subscript "Grhiup" is added to the coefficients named above 
with respect to equation 33 to indicate that the coefficients are for the high resistive state 
negative bias voltage. For equation 36, the subscript "Grhidn" is added to the coefficients 
named above with respect to equation 33 to indicate that the coefficients are for the high 
resistive state positive bias voltage. 

[0082] The values for the coefficients for equations 34, 35, and 36 for the data plotted in 
Figures 8, 9, and 10 are set forth in the column entitled "Values from Individual Curves" in 
Table 7 below. Also listed in the column are values for Ri 0 o and MR. The values in this 



-20- 



Attorney Docket No. SC13159TC 

column are utilized as coefficient parameter values in a model of an MRAM. The parameter 
values given in this column are fitted curve by curve. 

[0083] The last entry of this column, 105 ohms, is obtained, in one embodiment, by 
taking the root mean square error for the fit of the set of curves to the measured data. 



Table 7: Model Parameters 



Parameter 


Values from 


Globally 


Reduced 


Unit 




Individual 


Optimized 


Parameters 






Curves 


Values 






RloO 


16.9 


16.9 


16.9 


KOhms 


MR 


30.0 


30.5 


30.2 


% 


CjGrlo 


5.33E-4 


5.30E-4 


5.20E-4 


1/°C 


^VlGrlo 


5.62E-2 


5.81E-2 


6.21E-2 


1/V 


CrVlGrlo 


6.97E-5 


1.04E-4 


0 


i/v°c ! 


Cv2Grlo 


3.02E-1 


3.00E-1 


2.98E-1 


l/V 2 


CrV2Grlo 


5.35E-4 


6.26E-4 


6.57E-4 


l/V 2o C 


CjGrhi 


1.12E-3 


1.11E-3 


9.73E-4 


1/°C 


CviGrhiup 


-1.85E-1 


-1.91E-1 


-1.64E-1 


l/V 


CrviGrhiup 


3.51E-4 


6.79E-4 


0 


1/V°C 


Cv2Grhiup 


3.17E-1 


3.15E-1 


3.37E-1 


l/V 2 


CjV2GrhiuD 


7.61E-4 


1.02E-3 


3.93E-4 


l/V 2o C 


CviGrhidn 


3.49E-1 


3.58E-1 


3.51E-1 


l/V 


CyviGrhidn 


-4.53E-4 


-1.65E-4 


0 


irv°c 


Cv2Grhidn 


2.74E-1 


2.66E-1 


2.69E-1 


l/V 2 


CrV2Grhidn 


9.68E-4 


7.09E-4 


6.14E-4 


l/V 2o C 


Average Error 


105 


95 


96 


Ohms 



[0084] Referring back to Figure 12, in 1207, all of the parameters are combined to 
generate a global error measure for all of the curves. In one embodiment, the error measure is 
generated by the root mean square error, i.e. the square root of the average of the squares of 
the difference between the global model values and the individually measured actual 
conductance values. In another embodiment, the values at the different temperatures are 
given different weights before averaging. In one embodiment, the values nearer room 
temperature are given a higher weight than those farther away from room temperature. 



[0085] In 1209, the parameter values are adjusted to generate a least squares fit of the 
parameters to the combined data. The column entitled "Globally Optimized Values" of Table 
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7 sets forth the values of the parameters after then have been adjusted in 1209. This column 
shows that the root mean square error is reduced to 95 ohms 

[0086] In 1211, coefficient parameters which have a minimal affect on the global error 
are eliminated. In one embodiment, the parameters are eliminated if by setting them to zero 
and re-optimizing the other parameter values (e.g. performing operation 1209 again with the 
values shown in the column entitled "Reduced Parameters" in Table 7). In one embodiment, 
a parameter is eliminated (e.g. set to 0 in some embodiments) if the increase in global error is 
less than a particular threshold (e.g. %1) due to that parameter being eliminated. In some 
embodiments, re-optimization may be performed with a number of parameters set to zero. In 
other embodiments, the least effective parameters are set to zero until the cumulative effect 
exceeds a particular threshold. 

[0087] In some embodiments parameter candidates for re-optimization may be manually 
selected (e.g. by trial and error) from an evaluation of the data. In other embodiments, re- 
optimization is performed for with each parameter set to zero. In other embodiments, re- 
optimization is performed for each combination of parameters set to zero, wherein the case 
that eliminates the most parameters and yet provides an error below a threshold is selected. 

[0088] The column entitled "reduced parameters" of Table 7 shows a case where three of 
the parameters have been set to zero (GrviGrio, C T viGrhiu P , and OrviGrhidn) with the other 
parameters re-optimized. For this case, the total global error only increased by 1 ohm to 96 
ohms from 95 ohms for the globally optimized values prior to elimination. 

[0089] Providing a system for elimination of the parameters may provide for an MRAM 
cell model of increased efficiency in that less parameters are needed for calculation during a 
simulation of an MRAM memory. 

[0090] Figures 13-19 set forth one embodiment of a portion of a MTJ computer 
simulation model. The simulation model runs as a part of a simulation program for an 
MRAM memory which calls the model at each iteration (e.g. time step in some embodiments) 
of a simulation of an MRAM memory. In the embodiment shown, the flows set forth in 
Figure 13-19 are run for each time step (e.g. 100 picoseconds) of the simulation. 
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[0091] In the embodiment shown, the model implements the state machine of the 
Karnough Maps of Tables 1-6. The Karnough Map of Table 1 is implemented in the 
operations of Figure 13 from diamond 1309 on down, the Karnough Map of Table 2 is 
implemented in the operations of Figure 14, the Karnough Map of Table 3 is implemented in 
the operations of Figure 15, the Karnough Map of Table 4 is implemented in the operations of 
Figure 16, the Karnough Map of Table 5 is implemented in the operations of Figure 17, and 
the Karnough Map of Table 6 is implemented in the operations of Figure 18. 

[0092] Referring to Figure 13, there are two entry points for the model, initial entry 1301 
and other entry 1307. Initial entry 1301 is used when the model is called for the first time in a 
simulation in a simulation of a memory. Other entry 1307 is used on subsequent calls to the 
model including during subsequent time steps of the simulation. After initial entry 1301, 
various model parameters are initialized in 1303 including parameters for conductance as a 
function of the specific operating temperature and the initial bit state of the MTJ being 
modeled. In one embodiment, the model parameters are initialized (e.g. as set forth in Table 
7) by utilizing the globally optimized parameters or reduced parameters (if some parameters 
were able to be eliminated) (e.g. as set forth in Table 7) as derived by a process similar to that 
set forth in Figure 12. 

[0093] Once the initial conductance model parameters are set, the three conductance 
equations (low resistance state, high resistance state negative bias voltage, and high resistive 
state positive bias voltage) of an MTJ for simulation at a particular temperature looks like: 

G = G 0 [Cvo + Cvi V + C V2 V 2 ] (37) 

[0094] In one embodiment, the coefficient parameters C V o, Cvi, and C V 2 for each of the 
three conductance equations are calculated in 1303 using equations 34, 35, and 36, 
respectively, with the values derived by the method of Figure 12 as the coefficients for those 
equations and at the initial operating temperature. For example, Cvo (in equation 37 above) 
for the low resistance state case is calculated as (from equation 34) as being equal to l+C T Grio 
*(the difference between the reference temperature and the initialized temperature). 

[0095] In 1305, the model calculates the Axis 1 and Axis 2 currents (Ii and I 2 ) and the 
magnetic fields (Hi and H 2 ) for those fields from those currents. In the embodiment shown, 
Axis 1 is the write conductor located beneath the MTJ (e.g. metal line 1 14 in the embodiment 
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of Figure 1) and Axis 2 is the write conductor located above the MTJ (e.g. metal line 102 in 
the embodiment of Figure 1). 

[0096] The model then sequentially checks whether the magnetic field Hi generated by 
the current of Axis 1 (the Axis 1 magnetic field) has passed the lower threshold (Htl) while 
increasing (decision diamond 1309 of Figure 13) since the last iteration of the model, whether 
the Axis 1 magnetic field has passed the higher threshold (H T h) while increasing (decision 
diamond 1403 of Figure 14) since the last iteration of the model, whether the magnetic field 
H 2 generated by the current of Axis 2 (the Axis 2 magnetic field) has passed the lower 
threshold (H T l) while increasing (decision diamond 1503 of Figure 15) since the last iteration 
of the model, whether the Axis 2 magnetic field has passed the higher threshold (H T h) while 
increasing (decision diamond 1603 of Figure 16) since the last iteration of the model, whether 
the Axis 1 magnetic field has passed the lower threshold (Htl) while decreasing (decision 
diamond 1703 of Figure 17) since the last iteration of the model, and whether the Axis 2 
magnetic field has passed the lower threshold (Htl) while decreasing (decision diamond 1 803 
of Figure 18) since the last iteration of the model. 

[0097] If any of the checks is set forth in the preceding paragraph is positive, the model 
parses the state of the model according to the one Karnough Map in Tables 1-6 associated 
with the decision diamond that generated the positive check (e.g. the map of Table 1 is 
associated with diamond 1309, the map of Table 2 is associated with diamond 1403, the map 
of Table 3 is associated with diamond 1503, the map of Table 4 is associated with diamond 
1603, the map of Table 5 is associated with diamond 1703, and the map of Table 6 is 
associated with diamond 1803). 

[0098] Referring to Figure 13, if in diamond 1309, the model determines that the Axis 1 
magnetic field (Hi) has crossed the low threshold (Htl) while increasing, the model then 
determines in diamond 1311 whether a 2-1 pulse train has already started. The model 
determines whether the 2-1 pulse train has already started by determining whether the 2-1 
pulse train signal has been set (e.g. at a logical 1). If the 2-1 pulse train has been set, the 
model returns an error in operation 1313 in that such a condition is not logically reachable 
(e.g. it represents the state 0,0,1,0 of Table 1 which has an E). For the flows of Figures 13-19, 
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the "error exit" operations (e.g. 1313, 1317) represent states that should not be logically 
reachable. 

[0099] If no in diamond 131 1 (e.g. the 2-1 pulse train signal is a logical zero), the model 
then determines in diamond 1315 whether pulse 1 has already started. The model determines 
whether pulse 1 has started by determining whether the pulse 1 signal has been set (e.g. is at a 
logical 1). If yes in diamond 1315, the model returns an error in operation 1317 indicating an 
unreachable state. 

[0100] If no in diamond 1315, the model determines whether the 1-2 pulse train has 
already started in diamond 1327. The model determines whether the 1-2 pulse train has 
already started by determining whether the 1-2 pulse train signal has been set. If yes in 
diamond 1327, the model determines whether pulse 2 has already started in diamond 1323. 
The model determines whether pulse 2 has already started by determining whether the pulse 2 
signal has been set. If yes in diamond 1323, then the model sets the pulse 1 signal (to 1) in 
1321. If no in diamond 1323, the model returns an error in operation 1325. 

[0101] If no in diamond 1327, the model determines whether pulse 2 has already started 
in diamond 1333. If yes in 1333, the model determines whether the tau 1 pulse delay has 
been met in diamond 1331. If yes in diamond 1331, then the model sets the pulse 1 signal (to 
1), sets the 2-1 pulse train signal (to 1), and resets the tau 2 timer in operation 1335. If no in 
diamond 1333, then the model in diamond 1337 determines whether the tau 0 pulse delay has 
been met. If no in diamond 1337, then the model returns an error in operation 1339. If no in 
diamond 1331, the model returns an error in operation 1329. 

[0102] If no in diamond 1309, after operation 1321, after operation 1335, or if yes in 
diamond 1337, the model then goes to diamond 1403 of Figure 14 to determine whether the 
axis 1 magnetic field (Hi) has crossed the high threshold (Hth) while the magnitude has been 
increasing. 

[0103] If yes in diamond 1403, then the model determines in diamond 1405 whether the 
1-2 pulse train has already started and the 2-1 pulse train has already started. If yes in 
diamond 1405, then the model returns an error in operation 1407. If no in 1405, then the 
model determines in diamond 1409 whether pulse 1 has already started. If no in 1409, then in 
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diamond 141 1, the model determines whether the 1-2 pulse train has already started or the 2-1 
pulse train has already started in diamond 141 1. If yes in diamond 141 1, then the model 
returns an error in 1413. Ifnoinl411, then the model determines in diamond 1417 whether 
pulse 2 has already started. If yes in 1417, then the model returns an error in 1419. If no in 
1417, then the model determines whether the tau 0 pulse delay has been met in diamond 
1425. If no in diamond 1425, the model returns an error in operation 1427. If yes in diamond 
1425, the model sets the pulse 1 signal (to 1) and resets the tau 1 timer in operation 1431. 

[0104] If yes in diamond 1409, then the model in diamond 1415 determines whether 
pulse 2 has already started. If yes in diamond 1415, then the model determines in 1433 
whether the 1-2 pulse train has already started or the 2-1 pulse train has already started. If no 
in diamond 1433, then the model returns an error in operation 1435. If no in diamond 1415, 
then the model in diamond 1421 determines whether the 1-2 pulse train has already started. If 
yes in diamond 1421, then the model returns an error in operation 1429. If no in diamond 
1403, if no in diamond 1421, if no in diamond 1433, or after operation 1431, the model then 
goes to diamond 1503 of Figure 15. 

[0105] Referring to Figure 15, in diamond 1503, the model determines whether the axis 2 
magnetic field (H 2 ) has crossed the low threshold (Htl) while increasing. If yes in 1503, the 
model determines in diamond 1505 whether the 1-2 pulse train has already started. If yes in 
diamond 1505, the model returns an error in 1507. If no in diamond 1505, the model 
determines in diamond 1509 whether pulse 2 has already started. If yes in diamond 1509, 
then the model returns an error in operation 1511. If no in diamond 1509, the model 
determines in diamond 1517 whether the 2-1 pulse train has already started. If yes in 
diamond 1517, the model determines in diamond 1515 whether pulse 1 has already started. If 
no in diamond 1515, the model returns an exit in operation 1519. If yes in diamond 1515, the 
model sets the pulse 2 signal (to 1) in operation 1513. 

[0106] If no in diamond 1517, the model in diamond 1525 determines whether pulse 1 
has already started. If yes in diamond 1525, the model determines whether the tau 1 pulse 
delay has been met in diamond 1523. If no in diamond 1523, the model returns an error in 
operation 1521. If yes in diamond 1523, the model sets the pulse 2 signal (to 1), sets the 2-1 
pulse train signal (to 1), and resets the tau 2 timer in operation 1527. If no in diamond 1525, 
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the model in diamond 1529 determines whether the tau 0 pulse delay has been met. If no in 
diamond 1529, the model returns an error in operation 1531. 

[0107] If no in diamond 1503, after the pulse 2 signal has been set in 1513, after 
operation 1527, or if yes in diamond 1529, the model then proceeds to diamond 1603 of 
Figure 16. 

[0108] Referring to Figure 16, in diamond 1603 the model determines whether the axis 2 
magnetic field (H 2 ) has crossed the high threshold (Hth) while increasing. If yes in diamond 
1603, then the model in diamond 1605 determines whether the 1-2 pulse train has already 
started and the 2-1 pulse train has already started. If yes in diamond 1605, then the model 
returns an error in diamond 1607. If no in diamond 1605, then model determines in diamond 
1609 whether pulse 2 has already started. If no in diamond 1609, then the model determines 
whether the 1-2 pulse train has started or whether the 2-1 pulse train has started in diamond 
1611. If yes in diamond 1611, then the model returns an error in operation 1613. If no in 
diamond 1611, the model in diamond 1617 determines whether pulse 1 has started. If yes in 
diamond 1617, then the model returns an error in operation 1619. If no in diamond 1617, 
then the model in diamond 1625 determines whether the tau 0 pulse delay has been met. If no 
in diamond 1625, the model returns an error in operation 1627. If yes in diamond 1625, the 
model sets the pulse 2 signal (to 1) and resets the tau 1 timer in operation 1631. 

[0109] If yes in diamond 1609, the model determines whether pulse 1 has started in 
diamond 1615. If no in diamond 1615, the model determines whether the 2-1 pulse train has 
started in diamond 1621. If yes in diamond 1621, the model returns an error in operation 
1629. If yes in diamond 1615, the model in diamond 1633 determines whether the 1-2 pulse 
train has already started or the 2-1 pulse train has already started. If no in diamond 1633, then 
the model returns an error in operation 1635. If no in diamond 1603, if no in diamond 1621, 
after operation 1631, or if yes in diamond 1633, the model proceeds to diamond 1703 of 
Figure 17. 

[0110] Referring to Figure 17, in diamond 1703, the model determines whether the axis 1 
magnetic field (Hi) has crossed the low threshold (Htl) while decreasing. If yes in diamond 
1703, the model determines whether the 1-2 pulse train has already started and the 2-1 pulse 
train has already started in diamond 1705. If yes in diamond 1705, the model returns an error 
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in operation 1709. If no in diamond 1705, the model in diamond 1707 determines whether 
pulse 1 has already started. If no in diamond 1707, the model in diamond 1713 determines 
whether the 1-2 pulse train has started or the 2-1 pulse train has started. If no in diamond 
1713, then in diamond 171 1 the model determines whether pulse 2 has started. If yes in 
diamond 171 1, the model returns an error in operation 1715. If yes in diamond 1713, the 
model returns an error in operation 1717. 

[0111] If yes in diamond 1707, the model determines in diamond 1719 whether pulse 2 
has started. If no in diamond 1719, the model in diamond 1721 determines whether the 1-2 
pulse train has started. If yes in diamond 1721, the model returns an error in operation 1723. 
If no in diamond 1721, the model sets the pulse 1 signal and resets the tau 0 timer in 
operation 1731. After operation 1731, the model in diamond 1735 determines whether the 2- 
1 pulse train has started. If yes in diamond 1735, the model in diamond 1733 determines 
whether the tau 3 pulse delay has been met. If no in diamond 1733, the model returns an 
error in operation 1725. If yes in diamond 1733, the model resets the 2-1 pulse train signal in 
operation 1737. 

[0112] After operation 1737, the model in diamond 1759 determines whether the MTJ bit 
value has been set (whether the bit value is 1). If yes in diamond 1759, the model then resets 
the bit value (to 0) and sets the conductance parameters to the low resistive state values (e.g. 
using the coefficients for the low resistance state) in operation 1743. In one embodiment, the 
model sets the conductance parameters to low resistance state values to implement equation 
34 for determining the conductance provided by the MTJ in the circuit. 

[0113] If no in diamond 1759, the model sets the bit value of the MTJ cell (to 1) in 
operation 1747. After operation 1747, in diamond 1749, the model determines whether the 
bias voltage is positive. If in diamond 1749, the bias voltage is determined not to be positive 
(e.g. has a negative bias voltage) the model resets the down current flag and sets the 
conductance parameters to the high resistance state up values in operation 1755. In one 
embodiment, the model sets the conductance parameters to the high resistance state values 
such that it will use equation 35 in providing the conductance for the MTJ. 
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[0114] If yes in diamond 1749, the model sets the down current flag (to 1) to indicate that 
the bias voltage is positive in operation 1757. Also in operation 1757, in one embodiment, 
the model sets the conductance parameters to the high resistance state down values. 
Accordingly, in this embodiment, the model will use equation 36 when providing the 
conductance of the MTJ. 

[0115] If yes in diamond 1719, the model determines in diamond 1727 whether the 1-2 
pulse train has not started and the 2-1 pulse train has not started. If yes in diamond 1727, the 
model returns an error in operation 1729. If no in diamond 1727, the model in diamond 1739 
determines if the 1-2 pulse train has already started. If no in diamond 1739, the model in 
1741 resets the pulse 1 signal and resets the 2-1 pulse train signal. If yes in 1739, the model 
in 1745 determines whether the tau 2 pulsed delay time has been met. If no in diamond 1745, 
the model returns an error in operation 1751. If yes in diamond 1745, the model sets the 
pulse 1 signal and resets the tau 3 timer in operation 1753. 

[0116] If no in diamond 1703, if no in diamond 1711, if no in diamond 1735, after 
operation 1741, after operation 1743, after operation 1753, after operation 1755, or after 
operation 1757, the model goes to diamond 1803 of Figure 18. 

[0117] Referring to Figure 18, in diamond 1803, the model determines whether the axis 2 
magnetic field (H 2 ) has crossed the low threshold (H T l) while decreasing. If yes in diamond 
1803, the model in diamond 1805 determines whether the 1-2 pulse train has already started 
and the 2-1 pulse train has already started. If yes in diamond 1805, the model returns an error 
in operation 1809. If no in diamond 1805, the model determines in diamond 1807 whether 
pulse 2 has started. If no in diamond 1807, the model determines in diamond 1813 if the 1-2 
pulse train has started or the 2-1 pulse train has started. If yes in diamond 1813, the model 
returns an error in operation 1817. If no in diamond 1813, the model determines in diamond 
1811 whether the pulse 1 has started. If yes in diamond 1 8 1 1 , the model returns an error in 
operation 1815. 

[0118] If yes in diamond 1807, the model determines in diamond 1819 if pulse 1 has 
started. If no in diamond 1819, the model determines in diamond 1821 if the 2-1 pulse train 
has started. If yes in diamond 1819, the model returns an error in operation 1823. If no in 
diamond 1821, the model sets the pulse 2 signal (to 1) and resets the tau 0 timer in operation 
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1831. After operation 1831, the model in diamond 1835 determines whether the 1-2 pulse 
train has already started. If yes is 1835, then the model determines whether the tau 3 pulse 
delay has been met in diamond 1833. If no in diamond 1833, the model returns an error in 
operation 1825. 

[0119] If yes in diamond 1833, the model resets the 1-2 pulse train signal in operation 
1837. After operation 1837, the model in diamond 1859 determines whether the MTJ bit 
value is set (to 1). If yes in diamond 1859, the model resets the MTJ bit value (to 0) in 
operation 1843. Also in operation 1843, in one embodiment, the model sets the conductance 
parameters to the low resistive state values. Accordingly, in one embodiment, the model will 
use equation 34 for providing the conductance of the MTJ. 

[0120] If in 1859, the MTJ bit value is not set, then the model will set the MTJ bit value 
(to 1) in operation 1847. After operation 1847, the model determines in 1849 whether the 
bias voltage is positive. If the biased voltage is not positive in diamond 1849, the model 
resets the down current flag in operation 1855. Also in operation 1855, in one embodiment, 
the model sets the conductance parameters to the high resistive state up values. Accordingly 
in determining the conductance of the MTJ, the model will use equation 35. 

[0121] If the bias voltage is determined to be positive in diamond 1849, the model will set 
the down current flag (to 1) in operation 1857. Also in operation 1857, in one embodiment, 
the model will set the conductance parameters to the high resistive state down values. 
Accordingly, the model would provide the conductance of the MTJ as set forth in 
equation 36. 

[0122] If yes in 1819, the model then determines in 1827 if the 1-2 pulse train has not 
started and the 2-1 pulse train has not started. If yes in diamond 1827, the model returns an 
error in operation 1829. If no in diamond 1827, the model determines in 1839 if the 2-1 pulse 
train has already started. If no in diamond 1839, the model then resets the pulse2 signal (to 0) 
and resets the 2-1 pulse train signal (to 0) in operation 1841. If yes in diamond 1839, the 
model determines in diamond 1845 whether the tau 2 pulse delay has been met. If no in 
diamond 1845, the model returns an error in operation 1851. If yes in diamond 1845, the 
model sets the pulse 2 timer signal (to 1) and resets the tau 3 timer in operation 1853. 
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[0123] If no in diamond 1803, if no in diamond 181 1, if no in diamond 1835, after 
operation 1841, after operation 1853, after operation 1843, after operation 1855, or after 
operation 1857, the model then proceeds to diamond 1903 of Figure 19. 

[0124] Referring to Figure 19, once the switching threshold checks are completed 
(diamonds 1309, 1403, 1503, 1603, 1703, and 1803), the MTJ bit state is checked in diamond 
1903. This check is redundant if the bit has just been toggled (in operations 1843 or 1743) 
but is required if it has not been toggled. If the bit state is set (at 1) as determined in 1903, 
then the status of down current flag is checked in 1905 and if yes in 1905, the bias voltage is 
checked in 1907 to see if the bias voltage across the bit has changed to negative, necessitating 
a change in the bit conductance parameters to the high resistance state up values in operation 
1909. Also in operation 1909, the down current flag is set (to 1). 

[0125] If the down current flag is not set in 1905, then the model determines whether the 
bias voltage is positive in diamond 191 1, and if yes in 191 1, resets the down current flag (to 
zero) and sets the conductance parameters to the high resistive state down values in operation 
1913. 

[0126] If no in diamond 1903, if no in 1907, if no in diamond 1911, after operation 1909, 
or after operation 1913, the model performs operation 1915. In operation 1915, if there has 
been a change in conductance parameter values during the present iteration or previous 
iterations for a preset simulation time period (e.g. a 100 picoseconds of simulation time) the 
model does a transition between parameter values over time to prevent sudden changes in the 
conductance values for the simulation. 

[0127] Finally, the current through the MTJ is calculated from the voltage using the 
conductance parameter values and capacitance models (not shown in the Figures) in operation 
1917. The conductance parameter values utilized depends upon the most recent conductance 
values set and the transition simulation performed in 1915 if the parameters were changed 
during recent iterations. After the MTJ current is calculated, control is returned to other 
portions of the simulation. 
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[0128] In one embodiment, the flow set forth in Figures 13-19 may be implemented in 
VERILOG-A. In other embodiments, the model may have other configurations and/or other 
implementations. 

[0129] Although the Tables and flow charts described herein describe a model that checks 
for magnetic pulses crossing thresholds, a model made be based upon calculated current 
through the write lines crossing equivalent thresholds in other embodiments. With some of 
these parameters, the conductance coefficients would be scaled accordingly. 

[0130] Figure 20 is a block diagram of a computer 2020 in accordance with one 
embodiment of the present invention which may be used to execute the methods discussed 
herein. Computer 2020 includes a computer processor 2022 and memory 2024 operably 
coupled by a bus 2026. Memory 2024 may include relatively high speed machine readable 
media such as DRAM, SRAM, ROM, FLASH, EEPROM, bubble memory, etc. Also 
operably coupled to bus 2026 are secondary storage 2030, external storage 2032, output 
devices such as a monitor 2034, input devices such as a keyboard (with mouse) 2036, and 
printers 2038. Secondary storage 2030 may include machine readable media such as hard 
disk drives, magnetic drum, bubble memory, etc. External storage 2032 may include machine 
readable media such as floppy disks, removable hard drives, magnetic tap, CD-ROM, and 
even other computers, possibly coupled via a communications line. It should be appreciated 
that there may be overlap between some elements, such as between secondary storage 2030 
and external storage 2032. Executable versions of computer software 2033, such as, for 
example, software for performing the MRAM simulation described herein, can be written to, 
and later read from external storage 2032, loaded for execution into memory 2024, or stored 
on secondary storage 2030 prior to loading into memory 2024 and execution. Also, the 
MRAM cell simulation model may be stored in secondary storage 2030 or external storage 
2032. 

[0131] In the foregoing specification, the invention has been described with reference to 
specific embodiments. However, one of ordinary skill in the art appreciates that various 
modifications and changes can be made without departing from the scope of the present 
invention as set forth in the claims below. For example, any software taught herein may be 
embodied on one or more of computer hard disks, floppy disks, 3.5 M disks, computer storage 
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tapes, magnetic drums, static random access memory (SRAM) cells, dynamic random access 
memory (DRAM) cells, electrically erasable (EEPROM, EPROM, flash) cells, nonvolatile 
cells, ferroelectric or ferromagnetic memory, compact disks (CDs), laser disks, optical disks, 
and any like computer readable media. Also, the flow diagrams may also be arranged 
differently, include more or less operations, be arranged differently, or may have operations 
that can be separated into multiple operations that can be performed simultaneously with one 
another. Accordingly, the specification and figures are to be regarded in an illustrative rather 
than a restrictive sense, and all such modifications are intended to be included within the 
scope of present invention. 

[0132] In one embodiment, the invention includes a method for simulating a 
magnetoresistive memory device of a magnetoresistive random access memory (MRAM) 
having a first conductor, a second conductor, and a magnetic tunnel junction (MTJ). The first 
conductor is disposed substantially orthogonal to the second conductor. The MTJ is disposed 
between the first conductor and the second conductor. The method includes calculating a first 
current in the first conductor, calculating a second current in the second conductor, and 
detecting an indication of a transition of one of the first current and the second current across 
a threshold. The method also includes modifying a status of an operating condition of a 
plurality of operating conditions in response to the detecting the indication of the transition 
and outputting a bit state that is dependent upon a status of the plurality of operating 
conditions. 

[0133] In another embodiment, the invention includes a method of simulating a memory 
device of a magnetoresistive random access memory (MRAM). The memory device has a 
magnetic tunnel junction (MTJ) with multiple free magnetic layers. The method includes 
calculating an indication of a first magnetic field applied to the MTJ, calculating an indication 
of a second magnetic field applied to the MTJ, and detecting an indication of a transition of 
one of the first magnetic field and the second magnetic field across a threshold. The method 
also includes modifying a status of an operating condition of a plurality of operating 
conditions in response to the detecting the indication of a transition and providing an output 
bit state for the memory device. The output bit state is dependent upon a status of the 
plurality of operating conditions. 
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[0134] In another embodiment, the invention includes a method for simulating a 
magnetoresistive memory device in an integrated circuit magnetoresistive random access 
memory (MRAM) having a first conductor, a second conductor, and a magnetic tunnel 
junction (MTJ). The first conductor is disposed substantially orthogonal to the second 
conductor. The MTJ is disposed between the first conductor and the second conductor. The 
method includes calculating an indication of a first magnetic field applied to the MTJ. The 
first magnetic field is generated by current in the first conductor. The method also includes 
calculating an indication of a second magnetic field applied to the MTJ. The second 
magnetic field is generated by current in the second conductor. The method further includes 
detecting indications of transitions of the first magnetic field and the second magnetic field 
across one or more thresholds and providing a state machine having one or more state 
variables with transitions in the state machine being dependent upon detected indications of 
transitions of the first magnetic field and the second magnetic field and a state of the one or 
more state variables. 

[0135] In another embodiment, the invention includes a computer readable medium 
having stored instructions for simulating a magnetoresistive memory device of a 
magnetoresistive random access memory (MRAM) including a first conductor, a second 
conductor, and a magnetic tunnel junction (MTJ). The first conductor is disposed 
substantially orthogonal to the second conductor. The MTJ is disposed between the first 
conductor and the second conductor. The MTJ has multiple free magnetic layers. The 
computer readable medium includes instructions for calculating an indication of a first 
magnetic field applied to the MTJ. The first magnetic field being generated by current in the 
first conductor. The computer readable medium includes instructions for calculating an 
indication of a second magnetic field applied to the MTJ. The second magnetic field is 
generated by current in the second conductor. The computer readable medium includes 
instructions for detecting indications of transitions of the first magnetic field and the second 
magnetic field across one or more thresholds and instructions for modifying a status of 
operating conditions in response to detecting the indications of transitions of one of the first 
magnetic field and the second magnetic field across one or more thresholds. The computer 
readable medium also includes instructions for outputting a bit state that is dependent upon 
the status of the operating conditions. 
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[0136] While particular embodiments of the present invention have been shown and 
described, it will be recognized to those skilled in the art that, based upon the teachings 
herein, further changes and modifications may be made without departing from this invention 
and its broader aspects, and thus, the appended claims are to encompass within their scope all 
such changes and modifications as are within the true spirit and scope of this invention. 
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